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A volume-weighted measure for eternal inflation 
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I propose a new volume-weighted probability measure for cosmological "multiverse" scenarios 
involving eternal inflation. The "reheating-volume (RV) cutoff" calculates the distribution of ob- 
servable quantities on a portion of the reheating hypersurface that is conditioned to be finite. The 
RV measure is gauge-invariant, does not suffer from the "youngness paradox," and is independent of 
initial conditions at the beginning of inflation. In slow-roll inflationary models with a scalar inflaton, 
the RV-regulated probability distributions can be obtained by solving nonlinear diffusion equations. 
I discuss possible applications of the new measure to "landscape" scenarios with bubble nucleation. 
As an illustration, I compute the predictions of the RV measure in a simple toy landscape. 



I. INTRODUCTION 



In cosmological scenarios such as the "recycling uni- 
verse" [l| or the string-theoretic landscape 0, 0, 0] i the 
fundamental theory does not predict with certainty the 
values of "constants of nature," such as the effective cos- 
mological constant and particle masses. The cosmologi- 
cal observables may vary significantly between different 
causally disconnected regions of the spacetime. Hence 
one may only hope to obtain the probability distribution 
of the cosmological observables. Heuristically, one would 
like to compute probability distributions of the cosmo- 
logical parameters as measured by an observer randomly 
located in the spacetime. However, eternal inflation pro- 
duces an infinite volume in which possible observers may 
find themselves. Thus one runs into an immediate diffi- 
culty of defining a "randomly chosen" location within a 
noncompact space. 

Observers may appear only after reheating; the physics 
after reheating is tightly constrained by current experi- 
mental knowledge. The average number of observers pro- 
duced in any freshly-reheated spatial domain is a function 
of cosmological parameters in that domain. Calculating 
that function is, in principle, a well-defined astrophysical 
problem that does not involve any infinities. Here I focus 
on obtaining the probability distribution of cosmological 
observables at reheating. 

The set of all spacetime points where reheating takes 
place is a spacelike three-dimensional hypersurface [f| 0, 
0] called the "reheating surface." The hallmark feature of 
eternal inflation is that a finite, initially inflating spatial 
3- volume typically gives rise to a reheating surface having 
an infinite 3-volume, and even (potentially) to infinitely 
many causally disconnected pieces of the reheating sur- 
face, each having an infinite 3-volume (see Fig. [TJ. This 
feature of eternal inflation is at the root of several techni- 
cal and conceptual difficulties known collectively as the 
"measure problem" (see Refs. d i, M, G~l G3, for 
reviews and discussions of this problem). 

To visualize the measure problem, it is convenient to 
consider an initial inflating spacelike region S of hori- 
zon size (an "ii-region") and the portion R = R(S) of 
the reheating surface that corresponds to the comoving 




Figure 1: A 1+1-dimensional slice of the spacetime in an eter- 
nally inflating universe (numerical simulation in Ref. [3). 
Shades of different color represent different regions where re- 
heating took place. The reheating surface is the line separat- 
ing the white (inflating) domain and the shaded domains. 



future of S. If the 3-volume of R were finite, the volume- 
weighted average of any observable quantity Q at reheat- 
ing would be defined simply by averaging Q over R, 



(Q) = 



(1) 



where 7 is the induced metric on the 3-surface R. How- 
ever, in the presence of eternal inflation 1 the 3-volume 
of R is infinite with a nonzero probability X(<f>o), where 
4> = 4>o is the initial value of the inflaton field at S. The 
function X((f>o) can be computed in slow-roll inflation- 
ary models where typically X(4>o) I The ge- 
ometry and topology of the infinite reheating surface is 
quite complicated. For instance, the reheating surface 
contains infinitely many future-directed "spikes" around 
never-thermalizing comoving worldlines called "eternally 
inflating geodesies" Qj], EU • In a spacetime diagram 
such as Fig. [I] these spikes reach out to a timelike infinity 
since the eternally inflating geodesies never intersect the 
reheating surface. It is known that the set of spikes has 



1 Various equivalent conditions for the presence of eternal inflation 
were examined in more detail in Refs. QJ, O and 0. Here I 
adopt the condition that X(ij>) is nonzero for all <j> in the inflating 
range. 
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a well-defined fractal dimension that can be computed in 
the Fokker-Planck approach [l5| . 

For an infinite volume of R, the straightforward aver- 
age ([1} of a fluctuating quantity Q(x) over a noncom- 
pact reheating surface R is mathematically undefined. 
The average (Q) can be computed only after imposing a 
volume cutoff on the reheating surface in some way. A 
volume cutoff (or a "measure") is, in effect, a physically 
motivated prescription that makes volume averages (Q) 
well-defined. 

Volume cutoffs are usually implemented by restricting 
the consideration to a finite portion V of the reheating 
domain R. After imposing a cutoff, one computes the 
"regularized" distribution p(Q\V) of an observable Q by 
gathering statistics over a large but finite volume V. The 
final probability distribution p(Q) is then defined as 

p(Q) = lim p(Q\V), 

provided that the limit exists. A cutoff prescription is a 
specific choice of the compact subset V and of the way 
V approaches infinity when the cutoff is removed. It 
has been found early on (e.g. [(IHH) that p(Q) depends 
sensitively on the choice of the cutoff. Without a natu- 
ral mathematical definition of the measure, one judges a 
cutoff prescription viable if its predictions are not obvi- 
ously pathological. Possible pathologies include the de- 
pendence on choice of spacetime coordinates f20l l2ll , the 
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and the "Boltzmann 



"youngness paradox" [2 
brain" problem [2 

The presently viable cutoff proposals fall into two 
rough classes that may be designated as "worldline-based" 
and "volume-based" measures; a more fine-grained classi- 
fication of measure proposals can be found in Refs. [1,(3. 
Here I propose a new volume-based measure called the 
"reheating-volume (RV) cutoff." 



II. REHEATING-VOLUME CUTOFF 

In the RV cutoff, the reheating surface is not being 
restricted to an artificially chosen domain. Instead, one 
simply selects only those initial regions S that, by rare 
chance, evolve into compact reheating surfaces R having 
a finite, fixed volume Vol(-R) = V. The ensemble £y of 
such initial regions S is a nonempty subset of the ensem- 
ble £ of all initial regions S. The volume- weighted prob- 
ability distribution p(Q\£v) 01 a cosmological observable 
Q in the ensemble £y can be determined through ordi- 
nary sampling of the values of Q over the finite volume V. 
The RV cutoff defines the probability distribution p(Q) 
as the limit of p{Q\£y) at V — > oo, provided that the limit 
exists. 

To develop an approach for practical computations in 
the RV cutoff, let us first consider the probability density 
p(V; (f> )dV of having finite volume Vol(R) £ [V, V + dV] 
of the reheating surface R that results from a single H- 
region with initial value <p — <Po- This distribution is 



normalized to the probability of the event Vol(i?) < oo, 
namely 



p(V;^ Q )dV = Prob(Vol(i?) < oo) = 1-X((f> ). (2) 



The probability density p(V;<fio) is nonzero since 
X(4>o) < 1. I call this p(V;4>o) the "finitely produced 
reheated volume" (FPRV) distribution. This and related 
distributions constitute the mathematical basis of the RV 
cutoff. 

Below I will use the Fokker-Planck (or "diffusion") for- 
malism to derive equations from which the FPRV dis- 
tributions can be in principle computed for models of 
slow-roll inflation with a single scalar field. Generaliza- 
tions of p(V; <f>o) to multiple-field or non-slow-roll models 
are straightforward since the Fokker-Planck formalism is 
already developed in those contexts (e.g. [3ll. l32l]). 

Let us now define the FPRV distribution for some cos- 
mological observable Q at reheating. Consider the prob- 
ability density p(V, Vq r ; </>o, Qo), where tfio and Qo are 
values of <j> and Q in the initial 77-region, V is the total 
reheating volume, and Vq r is the portion of the reheat- 
ing volume where the observable Q has a particular value 
Qr. The distribution p(V, Vq r ; 4>q, Qo) as a function of 
Vq r at fixed and large V is sharply peaked around a mean 
value ( Vq r |y) corresponding to the average volume of re- 
gions with Q — Qr within the total reheated volume V. 
Hence, although the full distribution p(V, Vq r ; <po, Qo) 
could be in principle determined, it suffices to compute 
the mean value \Vq r \ v ). One can then expect that the 
limit 

P{QR) = J™ — ST" = v 1 ™ yProb(Vol(i^V) 

W 

exists and is independent of 4>q and Qo- (Below I will jus- 
tify this statement more formally.) The function p{Qr) 
is then interpreted as the mean fraction of the reheated 
volume where Q = Qr. In this way, the RV cutoff yields 
the volume-weighted distribution for any cosmological 
observable Q at reheating. 

To obtain a more visual picture of the RV cutoff, con- 
sider a large number of initially identical i7-regions hav- 
ing different evolution histories to the future. A small 
subset of these initial iJ-regions will generate finite re- 
heating surfaces. An even smaller subset of -regions 
will have the total reheated volume equal to a given value 
V. Conditioning on a finite value V of the reheating vol- 
ume, one obtains a well-defined statistical ensemble E\> 
of initial 7J-regions. For large V, the ensemble Ey can be 
pictured as a set of initial if-regions that happen to be lo- 
cated very close to some eternally inflating worldlines but 
do not actually contain any such worldlines (see Fig. 
In this way, the ensemble Ey samples the total reheat- 
ing surface near the "spikes" where an infinite reheated 
3- volume is generated from a finite initial 3-volume. It is 
precisely near these "spikes" that one would like to sample 
the distribution of observable quantities along the reheat- 
ing surface. Therefore, one expects that the ensemble Ey 
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Figure 2: A schematic representation of the ensemble E\> in 
comoving coordinates (t,x). Lightly shaded vertical strips 
represent the comoving future of various initial //-regions 
from E\>; dark shades represent reheated domains; the bound- 
ary of the dark-shaded domains is the reheating surface. Ver- 
tical dashed lines are the eternally inflating comoving world- 
lines that never cross the reheating surface. The 3-volumes 
of the reheating surfaces in the comoving future of the pic- 
tured //-regions are large but finite because these //-regions 
are located near eternal worldlines but do not contain any 
such worldlines. 



(in the limit of large V) provides a representative sam- 
ple of the infinite reheating surface, despite the small 
probability of the event Vol(R) = V. In this sense, the 
ensemble Ey at large V is designed to yield a controlled 
approximation to the infinite reheating surfaces R. 

The RV cutoff proposed here has several attractive 
features. By construction, the RV cutoff is coordinate- 
invariant; indeed, only the intrinsically defined 3- volume 
within the reheating surface R is used, rather than the 
3-volume within a coordinate-dependent 3-surface. The 
results of the RV cutoff are also independet of initial con- 
ditions. This independence is demonstrated more for- 
mally below and can be understood heuristically as fol- 
lows. The evolution of regions S conditioned on a large 
(but finite) value of Vol(i?) is dominated by trajectories 
that spend a long time in the high-// inflationary regime 
and thereby gain a large volume. These trajectories for- 
get about the conditions at their initial points and es- 
tablish a certain equilibrium distribution of values of Q 
on the reheating surface. Hence, one can expect that the 
distribution of observables within the reheating domain 
R will be independent of the initial conditions in S. 

The "youngness paradox" arises in some volume-based 
prescriptions because //-regions with delayed reheating 
are rewarded by an exponentially large additional volume 
expansion. However, the RV measure groups together the 
//-regions that produce equal final reheated volume V; a 
delay in reheating is not rewarded but suppressed by the 
small probability of a quantum fluctuation at the end of 
inflation. Therefore, most of these //-regions have "nor- 
mal" slow-roll evolution before reheating. For this reason, 
the youngness paradox is absent in the RV measure. A 
more explicit calculation confirming this conclusion will 
be given in Sec. IIVI 



III. RV CUTOFF IN SLOW-ROLL INFLATION 

As a first specific application, I implement the RV mea- 
sure in a slow-roll inflationary model with a scalar infla- 
ton d> and the action 



S 



16^ + \V^- V ® 



-gd 4 x. (4) 



If the energy density is dominated by the potential en- 
ergy V((j>), the Hubble expansion rate H is approximately 
given by 



H 



lirG 



V{cj>). 



(5) 



In the stochastic approach to inflation, 2 the semiclassical 
dynamics of the field <f> averaged over an //-region is a 
superposition of a deterministic motion with velocity 



<t> = v{<l>) 



1 



4ttG 



-M, 



(6) 



and a random walk with root-mean-squared step size 



(WY 



gW _ / ggw 

2tt V H{<j>) ' 



H 3 
Stt 2 



during time intervals St = H 1 . A convenient description 
of the evolution of the field at time scales St < H^ 1 is 



0(t + St) = <p(t) + v{cj>)8t + Z(t)^2D(cj>)St, (8) 
where £(t) is (approximately) a "white noise" variable, 

(o = o, mm) = s(t-t'), (9) 

which is statistically independent for different H- 
regions (36|. This stochastic process describes the evo- 
lution <j){t) along a single comoving worldline. For conve- 
nience, we assume that inflation ends in a given //-region 
when (j) reaches a fixed value <j> — 

Consider an ensemble of initial //-regions S\, S2, 
where the inflaton field <f> is homogeneous and has value 
4> = 4>q within the inflationary regime. Following the 
spacetime evolution of the field <j) in each of the regions Sj 
along comoving geodesies, we arrive at reheating surfaces 
Rj where <\> — 0*. Most of the surfaces Rj will have 
infinite 3-volume; however, some (perhaps small) subset 
of Sj will have finite Rj . The </>o-dependent probability, 
denoted X((f>o) = 1 — X((f> ), of having a finite volume of 
Rj is a solution of the gauge-invariant equation [l5j 



D{4>)X M + v(<t>)X t4 , + iH{4>)X lnX = 0, 



(10) 



2 See Refs. [H, [H, [H for early works and Refs. [TJ] and [TJj for 
pedagogical reviews. 
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with the boundary conditions = 1 and X(<pp{) = 1 

at reheating and at Planck boundaries. While X(<p) = 1 
is always a solution of Eq. (jlOp. the existence of a non- 
trivial solution with < X(<f)) < 1 indicates the pos- 
sibility of eternal inflation. The gauge invariance of 
Eq. (TlO|) is manifest since a change of the time variable, 
r(t) = J T(cj>)dt, results in dividing the three coefficients 
D,v,H by T{4>) which leaves Eq. (TO]) unchanged. 

The probability distribution p(V; 0o) can be found by 
considering a suitable generating function. Let us define 
the generating function g(z\ 4>q) by 



where D^, Dq, i^, vq are the suitable kinetic coefficients 
representing the "diffusion" and the mean "drift velocity" 
of <p and Q. The boundary condition at <fi — </>* is 

g(z, q; (j)*,Q) = exp [- (z + qS QQR ) iJ~ 3 (<?(>*)] , (16) 

where we use the delta-symbol defined by 5qq r = 1 if Q 
belongs to a narrow interval [Qr, Qb, + dQ] and 5qq r = 
otherwise. 

To obtain the distribution ([3]), we need to compute 
the average (Vq r | v ) at fixed V. We define the auxiliary 
generating function 



V<oo 



-zV 



p(V;0 o )dV. (11) h(z^,Q)^(V QR e- zV ) v<oo =-l q ( Z ,q = 0;^Q). (17) 



(Note that the formal parameter z has the dimension of 
inverse volume. The parameter z can be made dimen- 
sionless by a trivial rescaling which we omit.) The func- 
tion g(z; 4>o) is analytic in z and has no singularities for 
Rez > 0. Moments of the distribution p(V; (j>o) are de- 
termined as usual through derivatives of g(z; 4>o) in z at 
z = 0, while p(V, </>o) itself can be reconstructed through 
the inverse Laplace transform of g(z; </>o) in z. 

The generating function g(z;4>) has the following mul- 
tiplicative property: For two statistically independent 
-H-regions that have initial values <p — <pi and (f> = 4>2 
respectively, the sum of the (finitely produced) reheat- 
ing volumes Vi + Vi is distributed with the generating 
function 

(e- 2(Vl+V2) ) = ( e -^)<e-^)=.g(z;0 1 ). g (z;0 2 ). (12) 

This multiplicative property is the only assumption in 
the derivation of Eq. (flTJ]) in Ref. [15j. Hence, g(z;<p) 
satisfies the same equation (we drop the subscript in 



Dg,4>4> + v 9,<t> + 3Hg In g = 0. 



(13) 



The boundary condition at is g{z]<j) <t ) = e 
since an 7i-region with <p — 0» is already reheating and 
has volume i? -3 (</>*). The boundary condition at the 
Planck boundary (f>p\ (or other boundary where the ef- 
fective field theory breaks down) is "absorbing," i.e. we 
assume that regions with </> = <fip\ disappear and never 
generate any reheating volume: g(z;4>pi) = e z '° = 1. 
Note that the variable z enters Eq. ([13)) as a parameter 
and only through the boundary conditions. At z = 
the solution of Eq. (flU)) is g(0; 4>) = X(<p). Explicit ap- 
proximate solutions of Eq. (|13l) can be obtained using the 
methods developed in Ref. [la ]. 

Let us now consider FPRV distributions of cosmolog- 
ical parameters Q. The generating function for the dis- 
tribution piy, Vq r ; 4>, Q) discussed above is 

g(z, q;cj>, Q) = [[ e~ zV -^p(V, V Q ; 0, Q)dVdV Q . (14) 



The equation for g(z, q; <p, Q) is derived similarly to 
Eq. (fl~3|) and is of the form 



Note that g(z,q = 0; <fi, Q) = g(z]<p). The equation for 
h(z; (j), Q) then follows from Eq. (fT5|) , 

D^h t ^ + D Q h tQ Q +v 4> h^ + v Q h t Q + 3H (lng + 1) h = 0. 

(18) 

This linear equation contains as a coefficient the function 
g(z;(f>), which is the solution of Eq. (fT3"|) . The boundary 
condition for Eq. (fT8f is 



h(z- 0», Q) = e~ zH 3 ^H- 3 (<t>*)5(Q - Q R ). (19) 

Here we can use the ordinary ^-function instead of the 
symbol Sqq r because the 5-function enters linearly into 
the boundary condition. An appropriate rescaling of the 
distribution h by the factor dQ is implied when we pass 
from 6q Qr to 6(Q - Q R ). 

Finally, the expectation value {Vq r \ v ) at a fixed V 
and the limit ^ can be found using the inverse Laplace 
transform of h(z\ 0, Q) in z. 

The computation just outlined allows one, in principle, 
to obtain quantitative predictions from the RV measure. 
Further details and a more direct computational proce- 
dure will be given elsewhere [13] • Presently, let us ana- 
lyze the limit V — > oo in qualitative terms. The function 
p(Qr', V) is expressed as [cf. Eq. ©] 

(Vq b L) f°° e zV h(z;<P,Q)dz 
p(Q fl ;y)^^l^ J -'°? . (20) 



V 



Vjy^giz-^dz 



V$9M + D Q9,QQ + v 4>9,4> + V Q9,Q + ^Hglng = 0, (15) 



The asymptotic behavior of the inverse Laplace trans- 
form of h(z;<f),Q) at large V is determined by the loca- 
tions of the singularities of h(z; <fi, Q) in the complex z 
plane. The dominant asymptotics of the inverse Laplace 
transform are of the form oc exp(z*V), where z* is the 
singularity with the smallest |Rez*|. It can be shown 
that solutions of Eqs. (fT"3"|) and (fT5|) cannot diverge at fi- 
nite values of <p or Q. Thus g(z; 4>) and h(z; </>, Q) cannot 
have (f>- or Q-dependent singularities in z. Moreover, the 
function g(z; (f) cannot have pole-like singularities in z; 
the only possible singularities are branch points where 
the function g{z;<j>) is finite but a derivative with respect 
to z diverges. Furthermore, derivatives d™h satisfy lin- 
ear equations with coefficients depending on the deriva- 
tives d™~ 1 g(z; (p), which diverge at the singularities of g. 
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Hence the singularities of <3™/i in the z plane coincide with 
those of g(z; <j>). For these reasons the limit in Eq. ([3]) ex- 
ists and is independent of the initial values </>, Q. A more 
detailed analysis justifying these statements will be given 
in Ref. 



IV. RV MEASURE FOR A LANDSCAPE 

The "landscape" scenarios where transitions between 
metastable vacuum states occur via bubble nucleation, 
promise to explain the values of presently observed cos- 
mological parameters, such as the effective cosmological 
constant A. For this reason it is important to be able to 
apply the RV measure to landscape-type scenarios and, 
in particular, to compute the relative abundances of dif- 
ferent bubble types. 

In the terminology of Ref. [3|| , "terminal bubbles" are 
those with nonpositive value of A. No further transi- 
tions are possible from such bubbles because bubbles 
with A < rapidly collapse while bubbles with A = 
do not support tunneling instantons. The RV measure, 
as presently formulated, can be used directly for compar- 
ing the abundances of terminal bubbles. (Extending the 
RV prescription to non-terminal bubble types is certainly 
possible but is delegated to a future publication.) 

It was shown in Ref. [25| that the bubble volume cal- 
culations may use the simplifying "square bubble" ap- 
proximation, which neglects the effects of bubble wall 
geometry. In this approximation, the evolution of the 
landscape is well described by the "inflation in a box" 
model [lq . defined as follows. All the vacuum states are 
labeled by j — 1, N of which the terminal states are 
j = 1, Nt- During a time step St, an initial H- region 
of type j expands into rij = e 3HjSt independent daughter 
H- regions of type j . Each of the daughter H -regions then 
has probability Tj k of changing into an iJ-region of type 
k (for convenience we define Tjj = 1 — J2k^j ^jfe)- The 
process is repeated ad infinitum for each resulting H- 
region, except for //-regions of terminal types. A newly 
created //-region of terminal type will admit no further 
transitions and will not expand (or, perhaps, will expand 
only by a fixed amount of slow-roll inflation occurring im- 
mediately after nucleation). This imitates the behavior 
of anti-de Sitter or Minkowski vacua that do not admit 
further bubble nucleations. 

To implement the RV cutoff in this model, let us con- 
sider the probability Pj{n, n'; fc) of producing a finite to- 
tal number n of terminal //-regions of which n' are of 
type j, starting from one initial //-region of (nontermi- 
nal) type k. A generating function for this distribution 
can be defined by 



9j 



(z,q;k) = J2 z n q n ' Pj (n,n';k) = (z n q n ' 



n,n'—0 



following system of nonlinear algebraic equations, 



JV T N 

g 1 / nk (z : q;k)=Y,r kl zq^ + ]T T ki9j {z, q; *), (22) 



i=l 



I will merely sketch the derivation of Eq. (j2"2"|) , which is 
similar to the equations for generating functions used in 
the theory of branching processes (see e.g. the book [3^| 
for a mathematically rigorous presentation). The gener- 
ating function gj satisfies a multiplicative property anal- 
ogous to Eq. (TT2]) . This property applies to the n k inde- 
pendent daughter //-regions created by expansion from 
an //-region of type k. Therefore, gj(z,q;k), which is 
the expectation value of z n q n in an initial //-region of 
type fc, is equal to the product of n k expectation values 
of z n q n in the n k daughter //-regions (which may be of 
different types). The latter expectation value is given by 
the right-hand side of Eq. J32J. This yields Eq. after 
raising both sides to the power \jn k . 

If the generating functions gj are known, the distribu- 
tion pj(n,n';k) can be recovered by computing deriva- 
tives of gj{z, q; k) at z — and q — 0. Further, the mean 
fraction of //-regions of type j at fixed total number n of 
terminal //-regions is found as 



P(j\n) 



(n'\ n ) _ 0?0, gj (s = O,g=l;fc) _ 
n nd™gj(z = 0, q = 1; k) 



(23) 



Then the RV cutoff defines the probability of terminal 
type j, among all the possible terminal types, through 
the limit 



p(j) = lim p(j\n) 



+oo nd™gj(z = 0, q = 1; fc) 



lim 



(21) 

One can show that this generating function satisfies the 



(24) 

similarly to Eq. ([3]). Again one expects that the limit 
exists and is independent of the initial bubble type, as 
long as the initial bubble is not of terminal type. 

As a specific example requiring fewer calculations, let 
us consider a toy model with only three bubble types. 
There is one de Sitter (A > 0) vacuum labeled j = 3 that 
can decay into two possible anti-de Sitter terminal bub- 
bles labeled j = 1 and j = 2. The growth rate n 3 and the 
nucleation probabilities and are assumed known. 
To mimick interesting features of the landscape, let us 
also assume that there is a period of slow-roll inflation 
inside the bubbles 1 and 2, generating respectively N\ 
and N2 additional e-folds of inflationary expansion after 
nucleation. Hence, the model is determined by the pa- 
rameters 713, T31, T32, Ni, and A^. For convenience we 

define r 33 = 1 - r 31 - r 32 . 

We now perform the calculations for the RV cutoff in 
this simple model. There are only two generating func- 
tions, gi(z,q;k) and g2{z,q;k), that need to be consid- 
ered. The only meaningful initial value is fc = 3 (i.e., 
the initial bubble is of de Sitter type) since the two other 
bubble types do not lead to eternal inflation. Hence, we 
will suppress the argument fc in gj{z, q; fc). We also need 
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to modify Eq. to take into account the additional 
expansion inside the terminal bubbles. Let us denote the 
volume expansion factors by 

=3JVi 7.. 



(25) 



Then the functions g\ and g 2 are solutions of 

!/«3 _ y 31 z Zi " Zi ' r - — z 



9i 



+T 32 z" +r 33 5i, (26) 
r 3 iz Zl + T 32 z z *q z * + T 33 g 2 . (27) 



An explicit solution of these equations is impossible for a 
general n 3 (barring the special cases n 3 = 2, 3, 4). Never- 
theless, sufficient information about the limit (|24|1 can be 
obtained by the following method. Introduce the auxil- 
iary function F(x) as the solution F > of the algebraic 
equation 



F 33 F, 



(28) 



choosing the branch connected to the value F(Q) = 0. (It 
is straightforward to see that Eq. (|28|) . has at most two 
positive solutions, and that F(x) is always the smaller 
solution of the two.) The generating functions g\ and g 2 
are then expressed through F(x) as 

gi (z, q) = F(T 31 z Z iq Z i + T 32 z Z *), (29) 
g 2 (z, q) = F(r 31 z z > + Y 32 z z -q z -). (30) 



The limit (|24[) involves derivatives of these functions of 
very high order with respect to z. We note that the 
function F(x) is analytic; thus the functions g\ and g 2 
are also analytic in z. 

To evaluate the high-order derivatives, we need an el- 
ementary result from complex analysis. The asymptotic 
growth of high-order derivatives of an analytic function 
f(z) is determined by the location of the singularities of 
f(z) in the complex z plane. For instance, we may expect 
an expansion around the singularity z* nearest to z = 0, 
such as 



f(z) = c + ci (z- z*) s + 



(31) 



where s ^ 0,1,2,... is the power of the leading-order 
singularity, and the omitted terms are either higher pow- 
ers of z — z* or singularities at points z'^ located further 
away from z — 0. The singularity structure (|3ip yields 
the large-n asymptotics with the leading term 



dz n 



ci(— z*y 



r(n - s) 
T(-s) 



(32) 



This formula enables one to evaluate large-n limits such 
as Eq. fgjp . 

To proceed, we need to determine the location of the 
singularities of F(x). Since F(x) is obtained as an in- 
tersection of a curve F 1 /™ 3 and a straight line x + r 33 F, 
there will be a value x — x* where the straight line is tan- 
gent to the curve. At this value of x the function F(x) 
has a singularity of the type 



F(x) = F(x*) + Fiy/x — + 0(x — x*), 



(33) 



where F% is a constant that can be easily determined; 
we omit further details that will not be required below. 
The value of x* is found from the condition that dF/dx 
diverge at x — x*. The value of dF/dx at x ^ is found 
as the derivative of the inverse function, or by taking the 
derivative of Eq. 



dF 



1 



(34) 



dx 1 _r 33 

This expression diverges at the values 

F(x*) = (n 3 T 33 y^ , (35) 

(36) 



x* = (n 3 - l)T 33 F{x*) = T 33 ^ /(n ^_ iy 



Note that T 33 « 1 and n 3 ^> 1, hence is a constant of 
order 1. 

Rather than compute the limit (f2"3")l directly, we will 
perform an easier computation of the ratio of the mean 
number of bubbles of types 1 and 2 at fixed total number 
n of terminal bubbles, 



'(i) 



''(2) 



a?a g gi(z = o, g = i) 

d^d q g 2 (z = 0,q=iy 



(37) 



The derivatives d q g± and d q g 2 can be evaluated directly 
through Eqs. (|2"9"|) - (f3"0"]) . For instance, we compute d q g\ 
as 



9gi(z,q) 



dq 



= F'(T 31 z z ' + T 32 z z *)T 31 Z lZ z \ (38) 



9=1 



It is clear that the functions d q g\ and d q g 2 have a singu- 
larity at z = z* corresponding to the singularity x — x„ of 
the function F(x), where z„ is found from the condition 



r 32 ^f 2 



(39) 



Let us analyze this equation in order to esimate z*. If 
the nucleation rates T 3 i and T 3 2 differ by many orders 
of magnitude, we may expect that one of the terms in 
Eq. (|39|) . say T 3 iz Zl , dominates. Then the value z* is 
well approximated by 



r 3 i 



l/Z l 



(40) 



This approximation is justified if the first term in Eq. 
indeed dominates, which is the case if 



r 3 2 



1/Z 2 



£31 



1/Zi 



(41) 



If the reversed inequality holds, we can relabel the bubble 
types 1 and 2 and still use Eq. (|4*0|) . If neither Eq. (|4Tj) nor 
the reversed inequality hold, the approximation (|40[) for 
z* can be used only as an order-of-magnitude estimate. 
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To be specific, let us assume that the condition (j4"Tj) 
holds. 

We can now compute the ratio |37|) asymptotically for 
large n using Eqs. ([55]) and (f3"2"| . The singularities of the 
functions <?i and g 2 are directly due to the singularity of 
the function F. Then the dominant singularity structure 
of the function (1551) is found as 



dgi{z,q) 



dq 



i=i 2. -^(T 31 z^+T 32 z^) 



y/z — Z* 



(42) 

This fits Eq. (f3"Tj) , where f(z) = d q g\(z,q — 1) and s — 
— 1/2. After canceling the common n-dependent factors, 
we obtain 



p(2) 



lim 

n— >oo 



'■(1) 



(2) 



T 32 Z 2 z z i 



(43) 



This is the final probability ratio found by applying the 
RV cutoff to a toy model landscape containing two ter- 
minal vacua. Substituting the approximation (I40[) , which 
assumes the condition (|4ip. we can simplify Eq. (|4"3"]l to 
a more suggestive form 



p(i) r 31 Z! 

P(2) ~ 



1^32^2 



£31 

3 - ;); 



-I + Z2/Z1 



(44) 



We can interpret this as the ratio of nucleation proba- 
bilities r 3 i/r 32 times the ratio of volume expansion fac- 
tors, Z1/Z2, times a certain "correction" factor. As we 
have seen, the correction factor is actually a complicated 
function of all the parameters of the landscape. The cor- 
rection factor takes the simple form 



£31 



-1+Z2/Z1 



(45) 



only if the condition (|4"Tj) holds. 

We note that the result (j4"4"|) is similar to but does 
not exactly coincide with the results obtained in previ- 
ously studied volume-based measures. For comparison, 
the volume-based measures proposed in Refs. |38[ and 
[3 both yield 



"31^1 



P(l) _ 

p(2) r 32 z 2 ' 



(46) 



which is readily interpreted as the ratio of nucleation 
probabilities enhanced by the ratio of volume factors. 
The "holographic" measure [13], which is not a volume- 
based measure, gives the ratio 



pW 

P(2) 



£31 

r 3 2 



(47) 



that does not depend on the number of e-folds after nu- 
cleation. While the discrepancy between volume-based 
and worldline-based measures is to be expected, the cor- 
rection factor that distinguishes Eq. (|44|) from Eq. ([46]) 
is model-dependent and may be either negligible or sig- 
nificant depending on the particular model. 

As a specific example, consider bubbles that nucleate 
with equal probability, T 3 i = r 32 -C 1, but have very 
different expansion factors, Z\ 3> Z 2 . Then the condi- 
tion (|4ip holds and the "correction" factor is given by 
Eq. (|45[) , which is an exponentially large quantity of or- 
der Eg-, 1 . Qualitatively this means that the RV measure 
rewards bubbles with a larger slow-roll expansion factor 
even more than previous volume-based measures. 

On the other hand, if Z\ = Z 2 but r 3 i > T 32 , the 
"correction" factor disappears and we recover the result 
found in the other measures. 

To conclude, we note that the result (|4"3"| does not de- 
pend on the durations of time spent during slow-roll infla- 
tion inside the terminal bubbles, but only on the number 
of e-folds gained. This confirms that the RV measure 
does not suffer from the youngness paradox. 

So far we were able to apply of the RV measure to the 
comparison of the abundances of terminal vacua. One ex- 
pects that, with some more effort, the RV measure can be 
extended to arbitrary observables in landscape models. 
Further work will clarify the advantages and limitations 
of the new measure. 
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